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ABSTRACT 

We study the super-orbital modulation present in the Cygnus X-1 X-ray data, usually attributed to the pre- 
cession of the accretion disk and relativistic jets. We find a new, strong, 326±2 d period modulation starting 
in 2005, in Swift/BAT and RXTE/ ASM fight curves (LCs). We also investigate Ve/a5B/ASM and Ariel WASM 
arcfiival data and confirm the previously reported ^^290 d periodic modulation, and therefore confirming that 
the super-orbital period is not constant. Finally, we study RXTE/ ASM LC before 2005 and find that the pre- 
viously reported ^-^150 d period is most likely an artifact due to the use of a Fourier-power based analysis 
under the assumption that the modulation has a constant period along the whole data sample. Instead, we find 
strong indications of several discrete changes of the precession period, happening in coincidence with soft and 
failed state-transition episodes. We also find a hint of correlation between the period and the amplitude of 
the modulation. The detection of gamma-rays above 100 GeV with MAGIC in September 2006 happened in 
coincidence with a maximum of the super-orbital modulation. The next maximum will happen between 2 and 
14 of July 2008, when the observational conditions of Cygnus X-1 with ground-based Cherenkov telescopes, 
such as MAGIC and VERITAS, are optimal. 

Subject headings: binaries: general — X-rays: individual (Cygnus X-1) 



1. INTRODUCTION 



Cygnus X- 1 (iBowver et al.l 19651) is the best established can- 
didate for a stellar mass black-hole (BH). It is composed of 
a 21 ± 8 Mq BH turning around an 09.7 lab companion 
of 40 ± lOMfT, (Ziolkowski 2005) in a circular orbit of 5.6 
days (iBrocksopp et al. iil999a) . High resolution radio imag- 
ing has unveiled the pres ence of a highly colimate d, relativis- 
tic dStirUng et al.l l200ib . radiatively inefficient ( Gall o et all 
l2005h jet. The X-ray source displays soft and hard states 
and relatively frequent failed transitions between them. There 
are strong evidences of a high ener gy non-thermal compo- 
nent extending up to s oft gamma-rays JMcConnell et al.l2002l: 
ICadoUe Bel et al.ll200 6). The steady emission of gamma-rays 
above 100 GeV is strongly constrained by the observations 
with MAGIC which has obtained, however, a very strong ev- 
idence of an inten se, fast flaring episode at these energies 
(lAlbert et al. Il2007h . 

A ~5.6 d period modulation, attributed to the or- 
bital motion of the compact object around the compan- 
ion, has been obser ved at various wa v elengths by nu- 
merou s authors (e.g., Poolev et al. 1 11 9991: IBrocksopp et al. I 
" L998; LachowiczeLaLl |2Q0a). ~On 

a sup er-orbita l "^290 d period was 
claimed by iPriedhorskv. Terrel & Holt I (119831) on the soft 
X-ray data recorded by Ve/fl5B/ASM (1969-1979) and 
ArielV/ ASM (1974-1980). Later, a -150 d periodic vari- 
ability has been reported by various authors (iPoolev et al. I 
i999{'Brocksopp et al. "1999b^, 'Ozdemir & Demircan I |2001t 
Benlloch et al. 2001, 2004; Lachowicz et al. 2006) using dif- 
ferent data samples ranging between April 1991 and Novem- 
ber 2003. It must be noted, however, that other significant 
modulations w ith periods ~200 d and ^420 d have been 
also found (e.g. lBenlloch et al. 11200 U 120041: iLachowicz et al. I 
|2006). 
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the other hand. 



In this letter, we search the latest Cygnus X-1 X-ray data 
for periodic modulations. We also perform a critical revision 
of the previous results obtained from archival X-ray data. Fi- 
nally, we put our results in the context of a multiwavelength 
description of the source. 

2. DATA SAMPLES 

The data samples analyzed in this work are summarized in 
Table [T] They are available through the High Energy As- 
trophysics Archive Research Center (HEASARC). We use 
data from four different instruments, namely: Vela 5B/ ASM, 
Ariel y/ASM RXTE/ ASM and Swift/BAT. All data are aver- 
aged into one-day bins, except when explicitely stated. No 
period ic behavior is found in the X-ray data dur ing the soft 
state ( Wen et al. 1 119991: ILachowicz et al. 1 12006") and there- 
fore we analyze data corresponding to the hard state data 
only. The i nterval MJD 42338-42 829 is dominated by soft 
flare events (iLiang & Nolan |[l98 4l) and hence excluded from 
Ve/fl5B/ASM and ArielV/ ASM analyses. The soft state pe- 
riods during the operation of RXTE/ ASM are identified as 
those for which the ratio of count rates in the bands C (5.0- 
12.1 keV) and A (1.3-3.0 keV) is lower than 1.2 and the to- 
tal count rate exceeds the mean value by more than 4 stan- 
dard deviations. The mean and standard de viations are com- 
puted from the interval MJD 50660-50990 ("Lacho wicz et al. I 
2006). This excludes from the analysis the following periods 
(MJD): 50087-50327, 50645-50652, 51002-51026, 51369- 
51397, 51445-51625, 51776-51952, 52093-52584, 52762- 
52878, 52982-53092, 53198-53528, 53780-53872. After 
this, two long intervals dominated by hard state (samples A 
and B in Table [U are defined and studied separately. Based 
on the results for RXTE/ASM, the intervals 53414-53528 
and 53780-53872 are also removed from the analysis of 
Swift/BAT LC. 

3. ANALYSIS 

We search the different data samples for perio dic signals us- 
ing the Lomb- Scargle (L-S) test of uniformity (iLomb II 19761: 
IScargleiri982h . The chance probability is the probabihty of 
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TABLE 1 
Cygnus X-1 analyzed data 



Instrument/ 


Energy 


Operation 


Stai-t 


End 


Time 


Number of 


1^1 


'^m/2 


subsample 


[keV] 


time 


[MJD] 


[MJD] 


span [d] 


points (m) 


[d-M 


[d"l] 


Vela 5B/ ASM 


3-12 


May 1969-Jun 1979 


40368 


44042 


3675 


1097 


2.7x10"'' 


0.15 


Ariel WASM 


3-6 


Feb 1976-Feb 1980 


42830 


44292 


1464 


740 


6.8x10"* 


0.25 


RXTE/ASM 1 


2-10 


Sep 1996-Sep 1999 


50328 


51444 


1117 


913 


9.0x10"* 


0.41 


RXTE/ ASM 2 


2-10 


Jun 2005-May 2008 


53529 


54592 


1064 


909 


9.4x10"* 


0.43 


Swift/BAT 


15-150 


Jun 2005-May 2008 


53529 


54598 


1070 


829 


9.3x10"* 


0.39 



obtaining a certain L-S test value (zq) or larger out of a purely 
Gaussian noise sample, and is given by Ppre{z > zq) = e"^". 
When several frequencies are inspected, the post-trial proba- 
biUty, i.e. the probability to get a L-S test value zq or higher 
for at least one of the scanned frequencies, is given then by 
'Ppostiz > Zq) = 1 - [1 - Ppreiz > ^o)]" wherc n is the 
number of independent scanned frequencies. 

Given m data points, there is a discrete finite set of m/2 
independent frequencies. For the case of evenly spaced 
data there is a natural set of frequencies: Vk = ^ {k ~ 
1, . . . , m/2) where T is the time spanned by the data set. 
The values of the Fourier transform powers for the natural 
frequencies are independent of one another The data set does 
not contain enough information to search for periodicities be- 
low vi = 1/T or above i'mf2 = tti/2T. The time span, num- 
ber of data points and the maximum and minimum accessible 
frequencies for the different studied data samples are shown 
in Tabled 

For each investigated data sample we produce the peri- 
odogram, where — log]^Q(7'post) is represented as a function 
of the frequency. A prominent periodic component in the 
data is visible as a peak in the periodogram at the relevant 
frequency (I'peak)- We consider as significant those peaks for 
which the post-trial probability is lower than 10"^^, equiva- 
lent to a deviation of 5a from the Gaussian noise case. We 
scan all natural frequencies, with an oversampling factor of 
5. This means that we scan 5m/2 evenly spaced frequen- 
cies, from vi to v^/2- The oversampling does not increase 
the number of trials in the post-trial chance probability, since 
the number of independent frequencies remains constant, but 
increases the precision of v^^k- To estimate the error ( Aj/peak), 
we use the standard deviation of i/peak over 100 random d ata 
samples obtained by bootstrap (iDavison & Hinkleyll2006h of 
the original LC. Then, we fold the LC into a phaseogram us- 
ing the period 1/fpeak- The phaseogram is produced using 50 
bins, to ensure a smooth description of the waveform. The 
time of the phase (Tq) is determined from a fit to the LC 
using a Cosine function, where the value of the frequency is 
fixed to t'peak- In this way, Tq corresponds to the maximum 
of the fitting Cosine function (although not necessarily to the 
maximum of the waveform). The modulation amplitude (A) 
is defined as the ratio between the amplitude and the mean 
value of the Cosine function obtained from the fit. Finally, 
we remove the periodic component of frequency i^peak from 
the LC (a process called prewhitening). This is done by sub- 
tracting the deviations of the phaseogram from its mean value 
throughout the LC. 

We subsequently search for the next most prominent peak in 
the prewhitened LC, and follow the whole process described 
above in an iterative fashion. This process is stopped when 
the obtained fpeak has a post-trial chance probability larger 
then 10"6^ 



4. RESULTS 

We first search for periodic signals in ^wZ/f/BAT and 
RXTE/ ASM2 data samples, which correspond to the same 
epoch, and which are analyzed in this work for the first time. 
The results are shown in Table |2] The L-S periodograms for 
both LCs are shown in Figure [T] A very strong, dominant 
periodic signal with period 326±2 d is found in both data 
samples. The LCs are shown in Figure |2] The modulation 
is clearly seen by eye, which is reflected by the extremely low 
values of T'post. The previously reported ^^150 d modulation 
is not found in these data. The pre-trial chance probabilities 
for such a modulation are IQ-i-^ and IQ-^-^ for RXTE/ASM 
and Swift/BAil data, respec tively. Figure [T] sh o ws the ^ 150 d 
modulation as reported by iLachowicz et al. I (l2006l) overlaid 
with RXTE/ ASM2 LC, confirming that such a modulation 
does not describe well the data. The phaseograms correspond- 
ing to the 326 d period are shown in Figure [3] We see that 
hard and soft X-ray LCs are strongly correlated, with Pear- 
son's correlation factor r = 0.97. A second, also strong, 
component with period '-^lOOO d is present in both data sam- 
ples. This corresponds to a long-term modulation of the X-ray 
flux with respect to the 326 d oscillation, but cannot be estab- 
lished as periodic since the period is similar to the total time 
spanned by the observations. An alternative explanation to 
the ~1000 d period will be given below. Finally, a third com- 
ponent is seen in RXTE/ASM 2 data sample at period ^5.6 d, 
compatible with the orbital modulation, visible in the peri- 
odogram (Figure [Til even before prewhitening. This modula- 
tion is not seen in Swif t/BAT data, for which P pre. — 10"^. On 
a similar energy band, iPaciesas et al. I (Il997h claimed a mod- 
ulation compatible with the orbital period in CGRO/BAISE 
LC between April 1991 and September 1996 , but the method 
used lacks of a mathematical justification. .Brocksopp et al. I 
([l999b) did not find any evidence for the orbital modula- 
tion in CGRO/B ATSE data between May 1996 and Septem- 
ber 1998. Finally, Lachowiczet al. (2006), using the whole 
BATSE light curve, reported a deviation of ~ la from the 
Gaussian noise case, insufficient to establish the presence of 
the orbital modulation in the data. 

We have searched Vela 5B/ ASM and Ariel V/ ASM LCs 
for periodic modulations and found peaks at P=(276±3) d 
and P=(288±3) d, respectively (see Table [2]l, in agree- 
ment w ith the results obtained bylPriedhorsky. Terrel & Holt I 
d 19831) and ILachowicz et alTI (l2006l) . By comparison with the 
P=(326±3) d present in Swift/BATand RXTE/ASM 2 LC, this 
shows that the period of the super-orbital modulation is vari- 
able. The corresponding phaseograms are shown if Figure [3] 
(two lowermost panels). Even if the relative dispersion of the 
points is larger due to the lower sensitivity of Ariel V/ASM 
and Vela 5B/ ASM, the waveform is still visible. They fol- 
low a very similar shape as those found for 5w//f/BAT and 
RXTE/ ASM2, albeit for a different period, as one expects 
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TABLE 2 
Results for periodic signal search 



Sample 



^peak 

[d] 



To 
[MJD] 



A 
[%] 



7'pc 



Swi^/BAT 


326± 


2 


54027.3 


25 


10-87 




1030± 


50 


53707.2 


6 


10-10 


RX7'£/ASM2 


326± 


2 


54032.2 


29 


10-96 




990± 


30 


53670.7 


13 


10-41 




5.600± 


0.002 


53670.3 


5 


10-11 


RXr£/ASMla 


248± 


9 


50377.2 


29 


10-30 


RXTEI hSM\h 


123± 


3 


50761.1 


11 


10-10 


RXTE/ASM Ic 


168± 


4 


51177.6 


18 


10-24 


RXTE/ ASM 1 


5.602± 


0.002 


51117.1 


4.1 


10-^-1 


Ariel V/ASM 


276± 


3 


42865.8 


14 


10-16 


Vela 5B/ ASM 


288± 


3 


40187.4 


21 


10-6.8 



"E 0.35 ^ 




53600 53SO0 



Fig. 2.— Swift/BAT and RXTE/ASM2 LCs. The shaded area shows an 
interval of soft state, identified by the criteria exposed in Section |2] and 
not considered in the analysis. The thick, red curves are the fits by Cosine 
functions to each subsample (see Section |3). The thin, green curve repre- 
sents the ~ 150 d sup er-orbital modulation using the ephemeris reported by 
Lachowi cz et al. I <2006y The yeijical purple line marks the time of the TeV 
signal reported bv lAlbert et al. 1 120071) 




Fig . 1 . — Periodograms for Swift/BAT and RXTE/ASM 2 Cygnus X- 1 sam- 
ples, showing the post-trial chance probability as a function of the scanned 
frequency. The horizontal line marks the line corresponding to a post-trial 
probability of 10^^'^. The aiTow marks the orbital frequency. 



if the underlying physical process is the same. The corre- 
lation factor for Swift/BAT and Ariel V/ASM (Vela 5B/ ASM) 
phaseograms is r = 0.73 (r = 0.57). However, the modu- 
lation amplitudes are significantly lower than for the case of 
RXTE/ASM. This could have a physical explanation, but it 
could also happen if the periodic modulation was not present 
in part of the LCs, which can be certainly not excluded. On 
the other hand, we do not find evidence for the orbital period 
in Vela 5B/ ASM or Ariel V/ASM LCs"*. It is worth noting that, 
given Vela 5B/ ASM and Ariel V/ASM sensitivities, we do not 
expect to detect an orbital modulation with an amplitude of 
~5% as the one we see in RXTE/ ASM2 data. The mean rel- 
ative variance of the data points in the phaseogram (which is 
a good estimate of the measurement error) are 46% and 13% 
for Vela 5B/ ASM and Ariel V/ASM respectively. Both values 
are well above the 5% modulation which is hence difficult to 
detect. We have crosschecked this by analyzing Monte Carlo 
(MC) simulated LCs for Vela 5B/ ASM and Ariel V/ASM. We 
use the same sampling as the measured LCs and simulate 
a 5% amplitude modulation convolved with 46% and 13% 
point-to-point random fluctuations, respectively. The analy- 
sis of these LCs yields no significant peak. 

Finally, we have searched RXTE/ASM 1 data sample 
for periodic modulations, and found 5 significant peaks 
at P = (148±1), (188±2), (310±11), (475±11) and 
(5.598±0.003) d, all with chance probabilities lower than 
2Q-10 'pjjg laffgj- corresponds to the orbital modulation, 
whereas the other four seem to denote a complex power spec- 
trum. We stress that some of these peaks h ave been found 
in previous studies of the RX TE/ASM LC (Benlloc h et al. I 
1200 IL I2004t lLachowiczetal.li2006) . The understanding of 
the super-orbital modulation can be greatly simplified if we 

^^ We note that shorter integration times have been used for this search in 
the Ve/fl5B/ASM LC since, with 1-day bins, the minimum accessible period 
is P=6.7 d, as shown in Table[T] 
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Fig. 3.— From top to bottom: Swift/BAJ, RXTE/ ASM2, Ariel V/ASM 
and Vela 5B/ ASM phaseograms folded using period P and time values (P 
[d],To[MJD])= (326, 54027), (326, 54027), (276, 42866) and (288, 40187), 
respectively. The values of Tq are obtained from the fit of a Cosine function. 
Data points and error bars correspond, respectively, to the mean count rate 
and variance measured within each phase bin. The vertical, red, dashed line 
corresponds to the phase of the TeV signal reported by Albert et al. (2007). 




Fig. 4. — RXTE/ASMl Ught curve. The vertical shaded lines show the 
intervals of soft or failed transition states, identified by the criteria exposed in 
Section|2] which delimit the three subsamples (la, lb and Ic) used for further 
analysis. The red curves are the fits by Cosine functions to each subsample 
(see Section|3). 



consider that the period can change along the observation time 
in a discrete way. We have analyzed separately the data be- 
tween each two consecutive soft or failed transition states, 
i.e. three samples, namely la=[50328-50644], lb=[50653- 
51001] and lc=[51027-51368] (see Figure E]). We obtain a 
single significant peak in each of them (see Table |2]i. We 
prewhiten the leading frequency in each of the subsamples 
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and merge them to look for sub-leading frequencies. The only 
remaining significant peak corresponds to the orbital modula- 
tion (see Table|2li. Figure|4]shows RXTEI ASM 1 LC, overlaid 
with the result of the fits to independent Cosine functions to 
the three defined subsamples. The general agreement with 
the data is remarkably good. We have also generated three 
MC samples corresponding to the samplings of la, lb and Ic 
subsamples, and pure sinusoidal modulations with the periods 
found for each of them, convolved with the measured point- 
to-point fluctuations. Then we have merged the three samples 
together and analyzed the resulting LC. We obtain significant 
peaks at P = (147±1), (190±2) and (523±12) d, in surpris- 
ingly good agreement with the results of analyzing the real 
data. We note that this effect could be also responsible of the 
- 1000 d periodicity of the Swift/BAT and RXTE/ ASM 2 data, 
since they also contain a soft state episode that might have 
changed the period of the super-orbital modulation prior to 
MJD=53780. 

5. DISCUSSION 

We find that Cygnus X-1 displays a super-orbital modula- 
tion, with a period that changes, probably in a discrete way 
and in coincidence with soft or failed state-transition phases, 
over time scales ranging from a few hundred days to sev- 
eral years. According to our findings, the very much cited 
^150 d period is most probably an artifact of applying a 
(sometimes biased) Fourier-transform based analysis to a data 
sample where more than one consecutive period modulations 
are present. Since 2005, Cygnus X-1 shows a very powerful 
and stable super-orbital modulation with a period of 326 ± 2 d. 

The super-orbital modulation is usually at- 
tribute d to the precession of the accretion 
disk dPried horskv. Terrel & Holt ' '19 83|) a nd re lativistic 
jet (iRome ro. Kaufman Bernado & Mir abel I l2002h . as a 
result of the tid al forces ex erted by the companion star 
on a tilted disk dKatz I Il973h . A mechanism for keeping 
the disk tilted can be provided by radiation pressure warp- 
ing (|Pg tterson"1977'; 'Pri ngle lfT99a IWijers & Pringlelll999t 
lOgilvie & Dubus 2001). In the case of tidally forced preces- 
sion, the expected period Pprec depends on the outer radius 

Rq and inclinatio n of the disk 6 as Pprec oc Ro cos^^ 6 
(iLarwood I Il998h . Then, the longer the period the larger 
the precession angle, and hence also a larger modulation 
amplitude is expected. This is in agreement with our results 
for RXTE/ASM, where we have found four different super- 
orbital periods, which follow this tendency (see Table |2] and 



Figure |4|i. 

A different issue is why the precession of the disk pro- 
duces a modulati o n in t h e X-ray flux. Some authors (e.g., 
Lacho wicz et al. I 120061: llbragimov. Zdziarski & Poutanen I 
2007) have considered the possibility that the precession 
movement changes the optical thickness along the line of 
sight. They reject this possibility since it seems unlikely 
due to the fine tuning required to produce the observed 
modulation amplitude, which in addition should depend on 
the energy, which is not confirmed by our observations. The 
multiwavelenth data seem to support a scenario where the 
emission itself is anisotropic. The precession modulation 
is detected at similar times with identical periods in radio, 
soft and hard X-rays during the hard state. A unified picture, 
where the anisotropy is provided by the jet, has been proposed 
by Bro cksopp e t al. (1999b). The soft X-ray emission is 
produced in the disk via bremsstrahlung of thermal electrons, 
and are then up-scattered to higher energies via Compton 
scattering in the hot corona or at the base of the relativistic 
jet (which precesses with the disk). The acceleration of 
electrons along magnetic field lines in the jet produces the 
radio emission by synchrotron emission. During the soft 
state, the jet and corona dissappear and no modulation is 
observed. According to our findings, once the source goes 
back to the hard state, the reconstructed disk and jet have 
different kinematical properties, and the modulation period 
changes. It seems that failed transitions produce a similar 
effect. 

MAGIC detected a fast and i ntense episode of e mission 
of gamma-rays above 100 GeV (lAlbert et al. 1 12007^ during 
MJD=54003, albeit at the limit of the telescope's sensitivity. 
This happened in coincidence with the soft and hard X-ray 
maxima (Figures |2] and ^ and an unusually b right outburst 
detected with INTEGRAL (iMalzac et al. II2008I) . It is interest- 
ing to note that, according to the ephemeris shown in Table |2] 
the next passage for the precession maximum will happen at 
MJD=54655±2, i.e. between 6 and 10 of July 2008. The 
observational conditions of Cygnus X-1 with ground-based 
Cherenkov telescopes, such as MAGIC and VERITAS, will 
be optimal during those days. 
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